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Abstract 

We estimate how well we will know the parameters of solar neutrino oscillations after Kam- 
LAND and Borexino. The expected error on Am^ is few per-mille in the VO and QVO 
regions, few per-cent in the LMA region, and around 10% in the LOW region. The ex- 
pected error on sin^ 26 is around 5%. KamLAND and Borexino will tell unambiguously 
which specific new measurement, dedicated to pp solar neutrinos, is able to contribute to 
the determination of 6 and perhaps of Arn?. The present data suggest as more likely 
outcomes: no measurement, or the total pp rate, or its day/night variation. 
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1 Introduction 

The solar neutrino anomaly revealed in Homestake 
Kamiokande [D, Gallex [D and SAGE g has motivated 
the upgrades SuperKamiokande (SK) |l, GNO [§, and a 
new generation of experiments: SNO KamLAND ||] 
and Borexino H]. In the longer term, there are plans to 
attempt the real-time detection of pp neutrinos, thus cov- 
ering the whole solar neutrino spectrum ||l0| . Depending 
on the choice of experimental technique, it is hoped that 
future sub-MeV experiments will be able to measure some 
of the following properties of the solar neutrino flux at 
sub-MeV energies 

• total rate; 

• day/night variations; 

• "seasonal" variations; 

• energy spectrum; 

• total rate in neutral and charge-current (NC and CC) 
reactions. 

Today there are few disjoint best-fit solutions (usually na- 
med LMA, LOW, VO, . . . ) and these measurements could 

tOn leave from dipartimento di Fisica deU'Universita di Pisa and 
INFN. 



identify the true one. In fact, the survival probabilities 
Pee{Eu) for the present best-fit oscillations are not much 
different at E^, ~ 10 MeV where wc have more experimen- 
tal data, but are significantly different at lower E^. 

However, these sub-MeV experiments will presumably 
start after SNO, KamLAND and Borexino have already 
identified the true solar-neutrino solution and determined 
the solar-neutrino oscillation parameters. In this case, one 
should change the perspective and evaluate the potential of 
new experiments to improve on the measurement of solar 
oscillation parameters — to be contrasted with "to prove 
the occurrence of oscillations" . From this point of view, 
we answer the question in the title by determining how 
well near-future experiments are expected to contribute 
to these measurements. This fixes the minimal necessary 
accuracy of new sub-MeV experiments.|2] 

* We do not consider other possible reasons why one could be 
interested in sub-MeV solar neutrino experiments. They can be used 
to verify (or contradict) existing results. They could set bounds on 
exotic solutions of the solar anomaly (such as transitions into ster- 
ile neutrinos, into extra-dimensional neutrinos, into anti-neutrinos) 
or on neutrinos with exotic properties (such as a monster decay rate, 
or FCNC interactions, or magnetic moment, or else [p, l[ ) . They can 
check the existence of the MSW effect and the oscillation pattern 
fixed by more precise experiments, or to test solar model predictions, 
detecting possible short-scale time variations. However, helioseismol- 
ogy already provides accurate experimental information on the static 
properties of the sun. Finally, they could be used to demonstrate the 
validity of a new experimental method or technique. 
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The conclusions (see table 1 or fig. |l|) crucially depend 
on the true value of the oscillation parameters Am^ and 
9. A sub-MeV detector able to do different measurements 
but only with modest accuracy is never relevant to the 
measurement of oscillations parameters. Conversely, cer- 
tain specific sub-McV measurement could be relevant, if 
done precisely enough. Near-future experiments are able 
to cover fairly well all the possible cases (often measuring 
Am^ with great precision), and will indicate unambigu- 
ously which is the remaining relevant sub-MeV measure- 
ment. 



2 The near-future situation 

We focus mainly on the oscillations of three active neutri- 
nos, assuming that atmospheric oscillations do not affect 
solar neutrinos (in the standard notation, this amounts 
to have a large Am?j.„, and a small ^13) because the rel- 
evant experiments [0 | |[ | J[ |[ |l| 0, [l| [I6 
(except LSND 0, see page H) indicate that this is the 
physically relevant case. A three flavour analysis of solar 



(and atmospheric) data is no longer relevant [ p8| : reactor 
experiments directly demand a small 6*13 < 15° at 95% 
CL, so that solar oscillations are presently determined by 
the usual two parameters, Am^ and 9^ 

The continuous lines in flg. show the present global 
fit ||l9|, More or less acceptable fits can be obtained 
for a wide range of /S.m? and for a large mixing angle. The 
best-fit solutions have Xbcst ~ ^3 from 41 experimental in- 
puts and 2 free parameters; they lie in the LMA region. 



and have a Am around 10" 



The other solutions 



with smaller Am^ (named LOW, QVO, VO) are not sig- 
nificantly worse: they poorly fit the data where a 'solar 
anomaly' is present, i.e. the total rates, but satisfy well the 
bounds from data consistent with no oscillations, i.e. spec- 
tral distortions, seasonal and day/night variations. How- 
ever, the discrimination is not sharp: e.g., if GNO should 
decrease the central value of the Gallium counting rate 
down to 65 ± 5 SNU, 2a below the present value), the 
global fit would favor the LOW solution. 

In order to simplify the discussion it is useful to focus 
on few benchmark points that span the qualitatively dif- 
ferent, still allowed solutions. These points (denoted as 
point A. B, C, ...) are listed in table 1, and drawn as 'x' 
in fig. |l|. Some of them have a high Ax^ = — Xbost (i-^- 
are significantly disfavoured, roug hly at (Ax^)^/^ standard 
deviations: for example there is a ~ 7(t evidence for LMA 
versus no oscillations) . We retained them in order to have 
a conservative sampling of all possible cases. The points D, 
M, N, O have a significantly lower Ax^ in a solar-model- 
independent analysis, where the solar-model predictions 

t^ia could be large enough to give detectable effect in future solar 
neutrino experiments. But we will not consider this possibility in 
our analysis, because in this case S13 will be measured much more 
precisely by future long-baseline experiments. 



for the Boron and Beryllium fluxes are not used (see |£l|] 
for a precise description). 

Fig. |l| and table 1 illustrate the expected near-future 
achievements.|^| In the last column of table 1 we summa- 
rize which sub-MeV measurements could provide us with 
further, useful information. When a certain sub-MeV mea- 
surement is crucial (not very interesting), we mark it by 
a '!' ('?'). While Am^ can be often measured very accu- 
rately, the determinations of 9 may be less reliable. In- 
deed, the results on 9 depend strongly on solar-model pre- 
dictions; they could change e.g. if the Beryllium flux were 
lower than its predicted value. For points outside the LMA 
region, the error on 9 (and consequently on Pee) can be 
accurately estimated only after knowing the true results 
of near-future experiments, and could differ by a factor 2 
from the values quoted in table 1. 

The values of 6Pee in the penultimate column of ta- 
ble 1 are the 1 standard-deviation, near-future uncertainty 
on the survival probability of the total pp rate, as de- 
tected using ^ v'^e' scattering, with electron kinetic 
energy larger than Tg/ > 0.1 MeV (taking into account the 
smaller i/^^r NC cross section, the suppression in the total 
rate due to oscillations is roughly given by 0.8 Pee + 0.2). 
This quantity is important for sub-MeV experiments pO| , 
since it informs us on how accurately those based on elastic 
scattering (as Heron, Clean, Xmass, Genius, Hellaz, 
etc.) should measure the pp rate. Within the uncertain- 
ties, SPee is also relevant for experiments based on inverse 
/3-decay (as Lens, Moon, etc.). 

One should recall that existing Gallium experiments 
are sensitive to pp neutrinos (above 0.24 MeV), that would 
induce half of their neutrino events in absence of oscilla- 
tions. In the present perspective, we are lead to study how 
well the true rate of pp neutrinos can be reconstructed from 
the total Gallium rate, after subtracting the values of the 
other fluxes (Boron, Beryllium, . . . ) 



ptotal pothers 



as measured by the near-future (and present) experiments. 
Concerning the Boron flux, SK and SNO find that, within 
their accuracy, the survival probability of Boron neutri- 
nos is energy-independent. Therefore, we directly use the 
Boron rate measured at SNO with 8.5% error. The er- 
ror on the Boron contribution to the Gallium rate is ±1.5 
SNU, dominated by the 30% uncertainty on its Gallium 
cross section . 

Concerning the Beryllium and GNO fluxes, Borexino 
(and, maybe, KamLAND) should measure them. We as- 
sume that the error on Beryllium and GNO contributions 
to the total Gallium rate will again be dominated by the 

■fWhen the uncertainty is non Gaussian, one standard deviation 
errors have been replaced by one half of two standard deviation er- 
rors, if this gives a more conservative result. The 90% and 99% CL 
contours in our fits in fig. hi actually correspond to Ax^ = 4.6 and 
9.2, as obtained for two unknown parameters (Am^ a nd 9 ) using the 
Gaussian approximation, that is reasonably accurate [E0|. 
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Figure 1: The continuous lines show the present global fit of solar data in the (tan^ 0, Ato^/ eV^) plane. The coloured 
regions are the simulated near-future fits (after KamLAND and Borexino) for few benchmark points, marked as x and 
listed in table 1. The left (right) plot shows fits 90% (139% j confidence level. The true oscillation solution is identified 
in all cases (LMA, SMA, LOW, QVO, VO). Some simulated fits give very small regions: their main features can be 
read from table 1. 
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Table 1: A few benchmark points (ordered according to the value ofArn^), their present status, their near-future status, 
and the most relevant measurement left for sub-MeV experiments. 6 Pee is the error on the average survival probability 
of pp neutrinos. In absence of oscillations, near future experiments will indirectly measure Pee with ±0.06 error. 



~ 10% uncertainties on their Gallium cross section, so that 
the total error on will be ±2 SNU. We neglect the 

pep, hep and F contributions, that arc smaller than this 
error. The experimental error on RqI^^ is today ±5 SNU, 
and GNO should lower it down to ±(3 4) SNU. This 
means that the pp rate will be known within ±(10 15)% 
uncertainty, depending on its actual value (here, we as- 
sumed a large mixing angle). In this example, the average 
value of the survival probability for pp neutrinos (around 
0.5) will be known with a ±0.06 uncertainty. 

In the rest of the paper wc motivate and discuss in 
detail the results summarized in table 1 and fig. |^. 

LMA and EI 

The present data favour the LMA region, and do not sig- 
nificantly disfavour energy independent (EI) solar oscilla- 
tions, obtained for Am^ > 2 10"'' eV^ psl. This region 
should be fully covered by the KamLAND experiment 
which will detect reactor Pe using the reaction PeP e~^n. 
The accuracy of KamLAND has been discussed in |24, 
|l], m m and more importantly in 1^. As in |^ we 
assume (too conservatively?) that KamLAND will use a 
cut on visible energy Ee+ + rUe > 2.6 MeV (where Ee+ = 
Eii — mn+mp) in order to avoid the background due to am- 
bient Ve- We assume a rate of 550 events per year without 
oscillations, a 2% overall uncertainty on the reactor flux 
and no background. The regions in fig. surrounding the 
points B, C, D, E show how well KamLAND is expected 
to determine the oscillation parameters Am? and 9, af- 
ter three years of data-taking. The value of Am? can be 
measured accurately because the initial spectrum is well 
known, and the energy resolution is sufficient to show the 



first oscillation dip. 

However, if Am^ > 2 10"'' eV^ (the precise value de- 
pends on the energy resolution) KamLAND will only see 
averaged oscillations; thence, it will be unable to measure 
Am^ with good sensitivity ||2l|. Therefore, for point A 
we considered a new reactor experiment, with a baseline 
of 20 km, a rate of 3000 events/year above Ee+ + me > 
2.6 McV, and used three years of data for the estimate. 
This experiment has been named 'sub-KamLAND' in ta- 
ble 1 and in the following, because it is less demanding 
than KamLAND. | 

Despite the assumed 2% uncertainty on the initial Ve 
flux. Pee (and consequently the mixing angle 0) can be 
measured with an error less than 2%, z/an oscillation sig- 
nal is seen. In fact, the accuracy in the determination of 
Pee depends strongly on its value, being maximal when 
6 ^ 7r/4. It is possible to detect small deviations from 
maximal mixing. This can be understood in a qualitative 
way by considering an energy bin around the first oscil- 
lation dip, where Pee ~ 0. It will contain a number of 
events A'obs much smaller than the no oscillation predic- 
tion, A^o- Having assumed a negligible background, the 
survival probability in that bin can be measured with er- 
ror (5 Pee ^ VNohs/No- This feature is more pronounced 

§Here one crucial question is: if Am^ ^2 IQ-^eV^, do we need 
such a reactor experiment, or long-baseline experiments will do a 
better job? The answer is that even using a neutrino factory beam 
long-baseline cxperiments_jiall be unable to measure the solar pa- 
rameters accurately: see |2l[ for a discussion. The main reason is 
that Am^2 effects are entangled with effects due to 613. To disen- 
tangle them one need to compare some combination of observables, 
different from the one that gives the highest rate {ue ^ i'^ at a base- 
line of L ~ 700 km). Detailed studies of neutrino-factory capabilities 
assumed that the snls.r parameters will be accurately measured by 
other experiments p7| . 
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at the hypothetical sub-KamLAND than at KamLAND, 
where the are produced by various reactors with dif- 
ferent path- lengths, so that the first oscillation dip is par- 
tially averaged. (The situation would improve if Kam- 
LAND could reconstruct the direction of the neutrinos). 

We now discuss how sub-MeV experiments could im- 
prove on this situation. Solar neutrinos with energy 



2^/2GF^^e 



IMeV 



10-5 eV^ 



(1) 



(where N'^ is the electron density around the region of neu- 
trino production) do not experience the MSW resonance in 
the sun. Therefore, their oscillation probability is roughly 
given by averaged vacuum oscillations. Pee ~ 1 — ^ sin^ 26. 
In first approximation, this implies that sub-MeV exper- 
iments have nothing to tell about Am^, but could give 
information on 6. Note however that KamLAND will fix 
the value of Pee with an error less than ±0.02, as shown 
in table 1. Though the precision of sub-McV experiments 
is ultimately limited by the 1% solar model uncertainty 
on the pp flux, it seems unrealistic to aim at this level of 
accuracy, and even difficult to compete with KamLAND 
(or other upgraded reactor experiments). On the other 
hand, with a precise determination of Am^ and 6 from a 
reactor experiment, sub-MeV experiments could be used 
to finally test solar model predictions, as suggested long 
time ago |2^. 

Reactor experiments can only measure sin^(20): the 
discrimination between 9 and t: /2~ 9 has to be performed 
by relying on matter effects acting on solar neutrinos. The 
present data prefer 6 < 7r/4, with a few standard devia- 
tions significance (itsprecisc value depends on the actual 
value of 9 and Arn^ ) . Which new experiments are best 
suited for this issue? Neutrino fluxes are better known at 
lower E,j, but in the LMA region matter effects are larger 
at ~ 10 MeV than at E^ < MeV. At larger E^, the 
uncertainty on the Boron flux can be circumvented by the 
NC and CC measurements at SNO (and SK). The NC/CC 
ratio today prefers 9 < 7r/4 at ~ 1.5 standard deviations, 
and could provide a direct discrimination in the near fu- 
ture. If Am^ is in the lower part of the LMA region, earth 
matter effects at SK and SNO |^ could discriminate 9 
from tt/2 — 9. Matter corrections to the spectrum of Pe 
from a future supernova should also give a clear discrimi- 
nation: unlike v^, cross a resonance if > 7r/4. 

LOW 

The most promising LOW signals in near-future experi- 
ments are earth matter effects at Borexino and, maybe, 
at KamLAND In fact, after having excluded the 

LMA region, KamLAND could be converted into a solar 
neutrino experiment. In absence of oscillations, Borex- 
ino (KamLAND) is expected to have a signal rate of 59 



(280) events/day for the central value of the Beryllium, 
CNO, pep fluxes predicted by |^ and a background rate 
of 19 (130) events/day, known with 6b = 10% (50%j3 error. 
Even if the background were more uncertain, the seasonal 
variation of the solar neutrino flux (due to the excentricity 
of the earth orbit) would allow a measurement of the signal 
rate with an interesting accuracy |3^ . We concentrate 
on Borexino and consider only the signal due to the Beryl- 
lium line (46 events/day), after three years of running. 

At present, it is not clear if it will be possible to mea- 
sure the CNO and pep contributions, using a sufficiently 
background-free energy region above the Beryllium line. 
A pessimistic attitude would require to assume that it 
will be only possible to use those events with recoil e en- 
ergy between 0.25 and 0.8 MeV, generated by Beryllium 
(46 events/day), CNO (10 events/day), pep (2 events day) 
neutrinos |^ plus of course the background. However, we 
remark that this possible limitation would be a problem 
to test of solar models, but would have just a little effect 
on the determination of Arrv^ and 9. 

We perform our analysis along the lines of js^, but 
paying more attention to the determination of the oscil- 
lation parameters, rather than to the discovery of oscilla- 
tion signals. We divide the simulated data into A'^bins = 
8-1-1 zenith bins (8 night zenith bins equally spaced in 
cos6'zenith, plus ouc day bin) and define 



^'Bc ^^Bc 
' Be 



(2) 



where Ni is the measured number of events; Si is the ex- 
pected number of signal events (proportional to Pee and to 
the solar flux ^7^^); hBi is the expected number of back- 
ground events (the factor h takes into account the overall 
uncertainty on the background rate). Finally the in 
eq. (H) is added to the from present data, properly tak- 
ing into account the correlation of the theoretical uncer- 
tainties on the Beryllium flux ||3l|, |2^. A ^s in eq. 
when minimized with respect to the 'nuisance' unknown 
parameters (here b and <&7Bg; more generically the solar- 
model parameters and the detection cross sections) gives 
the same defined in p^ , where all uncertainties are 
summed in quadrature, obtaining a big error matrix. 

The simulated fit for point G (which gives the best-fit 
in the LOW region) is similar to the corresponding result 
in fig. 6 of |32] (where KamLAND instead of Borexino 
was considered). The accuracy is worse at point H, be- 
cause earth matter effects diminish with Am? . We do not 
show simulated fits for points located around the high- 
est Am? values allowed in the LOW region (that will be 

^Here we hope to be too pessimistic. We assume that the fiducial 
volume will be 60% of Kami, AN D M j ; this is why our rates are 60% 
lower than those employed in [Bu, 32| . 
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soon tested by SNO and are already directly disfavoured by 
the non-observation of a day/night asymmetry at SK and 
Gallex/SAGE/GNO [||): due to the large earth matter 
effects, the accuracy on the determination of Am^ would 
be so good that one should carefully take into account the 
uncertainty in the profile density of the earth (while this 
is not an issue for the points that we have selected). 

The survival probability of sub-MeV neutrinos is given 
by adiabatic conversion: P^e = sin^ 9 during the day. An 
accurate measurement of the pp rate would provide the 
safest determination of 0, because solar-model predictions 
will play little role. It will be interesting to perform this 
measurement even if present and near-future experiments 
will nominally give a somewhat more accurate determina- 
tion of 9. 

Furthermore, earth matter effects give a day/night vari- 
ation of the pp rate, allowing to measure also Am? . How- 
ever, even if matter effects are larger at pp energies than 
at higher energy, it is more convenient to study Am^, 9 at 
Borexino (KamLAND, or new experiments based on in- 
verse /3-decay) due to the monochromaticity of Beryllium 
neutrinos and to the larger event rate. 

LOW/QVO boundary 

We pragmatically define the boundary between the LOW 
and QVO regions as 

Am2(L0W) > 10"^ eV^ > Am^iQYO) 

because this is the critical Am^ under which Borexino 
(KamLAND) should observe anomalous seasonal effects 
(see e.g. ||, |3^), rather than the day/night effects 
characteristic of the LOW region (see e.g. I^, ^). The 
regions considered in the rest of this paper will be soon 
disfavoured, if SNO finds a day/night asymmetry. 

No unmistakable signal of solar neutrino oscillations 
can be observed by near-future experiments if Am^ lies 
around this critical value. In view of this situation, we 
tried to exploit, in this particular point, the full capabil- 
ity of near-future experiments by optimistically assuming 
that KamLAND will be converted into a solar neutrino ex- 
periment as described above. KamLAND would detect a 
hint of day/night effect, giving some information on Am^. 
Borexino (and existing data) would provide instead the 
dominant information on 9. By combining these pieces of 
information, we obtain the simulated fit for point I. 

Though the near-future uncertainty in Am^ is signif- 
icantly smaller than the present uncertainty, it remains 
rather big. It should be understood that the improvement 
is mainly due to the fact that all other regions with larger 
and smaller Am? will be firmly excluded, because they 
predict unobserved clear signals. To prove this, we omit- 
ted the only positive signal (attributed to KamLAND) and 
obtained roughly the same Am^ interval. 



A measurement of the pp rate could give additional 
informations on 9. However a measurement of Am^ re- 
sulting from a global fit may be felt as unsatisfactory. As 
discussed in the 'no oscillations' section, the detection of 
Beryllium neutrinos would ensure that a solar neutrino 
anomaly is present, but we would still not know if it is 
due to oscillations. A sub-MeV experiment could discover 
an unmistakable oscillation signal: the neutrino energy at 
which earth matter effects induce a maximal day/night 
asymmetry is 



2V2GFiV| 



40 10-8 eV^^ 



where Nf is the electron density of the earth mantle. Be- 
cause of the low value of E^^^, one would need a big real- 
time detector, with an energy threshold as low as possible 
in order to detect this effect. 

However, at a detector with a threshold Tg > 0.1 MeV 
(such as Heron, Clean, Xmass, Genius, Hellaz etc.) 
the day/night asymmetry at point I (a ^ 4 % excess in 
nighttime) is only twice larger than at Borexino or Kam- 
LAND, that could instead have many more events. A very 
low threshold, maybe as low as 11 keV, could be attained 
by the Genius experiment with a rate of 18 pp events/day 
(assuming a mass of 10 ton). However only few % of the pp 
neutrinos have energy below 0.1 MeV, and they cannot be 
individually identified by the elastic scattering reaction. 



QVO 

The distance between the earth and the sun varies as 

L{t) = Lo[l - ecos(27rt/yr) + 0{e^)] 

where = 1-496 10* km is the astronomical unit, e = 
0.0167 is the excentricity of the earth orbit, and t is the 
time since the perihelion (around 4th of January). This 
variation induces a modulation of the survival probability 
in QVO and VO regions, that can be investigated at Borex- 
ino, and eventually at KamLAND, by means of the almost 
mono-energetic Beryllium neutrinos, Et^^ = 0.863 MeV [p5[ 



The time variation of the survival probability is |34 

Pee(i) = {Pee) + sin2(?V^c(l - Pc)D COs{kL{t) + 5) 



where k = Am^ jlEi^^ and (5 - O.13/fci?0 Q. The num- 
ber of oscillations met in one semester can be large (fig. ^): 



tLnAvr? 



27r£;7Bc 0.85 10-9 eV^ 



(3) 



These oscillations get washed for increasing Arv? , due to 
the MSW effects inside the sun (taken into account by the 
crossing-probability factors Pc) and to the finite width of 
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Figure 2: Number N^sc o-'<^d amplitude a of seasonal 



oscillations as function of Am 



{l/vTO, 1, a/IO} (long dashed, dotted, dashed). 



eV^ for tan^e 



Figure 3: Seasonal variation of the signal at Borexino for 
the benchmark points. The upper line is the no-oscillation 



the Be line (taken into account by the factor D, given 
by the Fourier transform |Q of the BeryUium 'line' spec- 
trum |3^), as illustrated in fig. 

We perform the simulated fits using a with a large 
number of seasonal bins: A'bins ~ 10 • iVosc- Of course, 
A'osc is a priori unknown in actual analyses, and should 
be extracted from the data; a definition of that avoids 
this problem, and allows to exploit all the data is discussed 
in Appendix A. Anyhow, our conclusion is that Am^ can 
be measured with surprisingly good precision (sec table 1, 
or enlarge fig. |^). The point is that when the number of 
oscillations A'osc is big, Borexino acts as an interferometer. 

These results can be understood by a simple analytical 
estimate. First, we display in fig. |^ the signal as a function 
of time, for the benchmark points I. L and M. The number 
of events at Borexino can be written as: 



A^ cx 1 + acos{kL + S). 



(4) 



considering the events due to Beryllium neutrinos (CC and 
NC) and the estimated background, and omitting the ge- 
ometrical flux factor. The amplitude of oscillations 
a is plotted in fig. ^, for three choices of the mixing an- 
gles. The average value of A^ during the time periods with 
cos{kL + (5) > (< 0) is N± cx 1 ± 2a/7r. Thence, one gets 
an "asymmetry" (systematic excess) of events: 



4a 

TT 



A+ - A^_ 



± 



(A'+ + A_)/2 yA+ + A_' 
The error on a can be small even if A'osc is big: 



6a 



N+ + N_ 



0.005, 



/ 50000 

A+ + A^_ 



The 'seasonal' oscillation is detected if the amplitude a is 
sufficiently larger than its uncertainty Sa, and this hap- 



pens if Ato2 < 0.5 10"** eV^ Values of 61 < 7r/4 are more 
suitable; sec again fig. |^. 

Once a seasonal signal is seen with A"o- = a /6a standard- 
deviations, the number of oscillations A'osc can be mea- 
sured with an error (5 A'osc ~ ^/'2Ncr that does not depend 
on Aosc- Using eq. (^), the consequent error on Am^ is 



SAm^ 
Aw? 



5N„. 



1 



A^o, 



if A^ > few. (5) 



For example, from fig. ^jwe can read that in point L Aqsc ~ 
7 and N„ ~ 20. This explains why in table 1 we claim an 
extremely accurate measurement of Ann?. The analytical 
approximation tells how our numerical results should be 
rescaled in order to extend our analysis to KamLAND or 
to a longer data-taking period. For example, with 9 times 
more signal, 5a would be 3 times smaller, giving a per- 
mille determination of Am^. The accuracy is limited by 
statistics (and by the small cxccntricity of the earth orbit). 
No knowledge of the Beryllium spectrum or of other solar- 
model dependent features is required. 

The overall phase in cos(A:io + S) depends on the total 
earth-sun distance (rather than on its excentricity varia- 
tion) and is therefore more strongly dependent on Am?. 
It contains additional information: it separates the al- 
lowed region of Am^ into (5 A'osc /2e ~ few separate thin 
islands |Q. We do not include this information in our 
fits: the error or Am^ is already so small that it would not 
even be possible to see this sub-structure in fig. |l|. 

Beryllium neutrinos contribution is ~ 30% of the total 
rate measured at SAGE and Gallex/GNO. These exper- 
iments are further limited by a modest time resolution. 
At ^ weeks, and thence their sensitivity to seasonal vari- 
ations is reduced. No such a signal has been found, which 
implies weak bounds on the oscillation parameters ||37|] , 
usually neglected in global analyses. 
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Subsequent pp experiments cannot improve on this de- 
termination of Am^. A measurement of the pp rate would 
instead give a useful information on the mixing angle. 
However, this information would have a certain degree of 
solar-model dependence. Indeed, in the QVO region the 
survival probability (Pee) = | + {Pc ^ 5) cos2^^ lies some- 
where between vacuum oscillations {(Pee) = 1 — ^ sin^ 20 
for Pc = cos^ 6) and adiabatic oscillations {{Pee) = sin^ 9 
for Pc = 0), as controlled by the crossing-probability Pc = 
[e^cos^ 9 _iy[g'r _ 1] .^^^ere Q 

^ ~ E^ldlnNe/drl-cs ^ lO-^eVVMcV' 

The gradient is evaluated around the resonance point (for 
a more accurate approximation see [ p4[ ) where the density 
is Ne ~ Am'^ /GpEi,: this corresponds to the outer part of 
the sun where the profile density deviates from the simple 
exponential approximation, A^e c>c exp(— 10.54r/i?sun)- 

Before concluding, we recall that an accurate treatment 
of CNO and pep neutrinos would require to know the ac- 
tual performance of the near-future detectors. However, 
the sensitivity of the determination of Am^ through Beryl- 
lium neutrinos would be not substantially affected, even in 
the most pessimistic case. Instead, in the most optimistic 
case, pep neutrinos could give an additional modulated sig- 
nal with frequency fcpep ~ 0.6k and amplitude Opep ^ a/10, 
referred to those of Beryllium neutrinos, eq. (^. Since the 
pep flux is well known, almost as the pp flux, a^ep could 
yield information on 9. 

vo 

Vacuum oscillations offer the same seasonal signal as dis- 
cussed in the previous section. Solar matter effects can now 
be neglected, so that the energy-averaged survival proba- 
bility is {Pee) = 1^5 sin^ 29. The number of seasonal 
vacuum oscillations A'osc encountered during one semester 
is smaller, so that the accuracy in Am^ is somewhat worse 
(see eq. (||)), but remains much better than sub-MeV ca- 
pabilities. A measurement of the pp rate would give a 
solar-model independent information on the mixing angle. 
Some vacuum oscillation 'solutions', usually named "Just 
So^" , that give poor fits of existing dataQ present charac- 
teristic spectral distortions in pp neutrinos . Our point 
N gives a rather good fit to the data. In fact, the Beryl- 
lium line is affected by seasonal effects (while the other 
fluxes are less affected due to their larger energy spread) 
in a way that very strongly depends on Am^: choosing 
the appropriate value one can obtain the experimentally 
preferred value of the Beryllium rate. 

*We get Ax^ ~ 15, which is in agreement with the analyses in | |l9| , 
except the second one. In these 'solutions', Pee deviates from unity 
only for low energy neutrinos: if the 3 standard-deviation evidence 
from SK and SNO against this possibility is correct, new SNO data 
will soon make it even more disfavoured. 



SMA 

The SMA region is strongly disfavoured by existing experi- 
ments and has a low goodness-of-fit. The reason is that the 
SMA oscillations that fit the measured rates imply a sur- 
vival probability Pee{E^) in conflict with the SK spectral 
data. This conflict can be seen in pre-SNO flts performed 
by the SK collaboration ||^ and has become sharper after 
SNO. Our SMA point, F, is not the 'best' current SMA 
solution (that has Ax^ « 14, simflar to ||l|, |o)) but a rep- 
resentative point of this region. In view of this situation, 
it is difficult to study seriously how well new experiments 
can measure Arri^ and 9 in the SMA region. 

The most characteristic feature of the SMA region is 
that 

1. the neutrino rate at Borexino will be strongly sup- 
pressed (almost down to the background level). 

Depending on the actual SMA oscillation, this evidence for 
SMA can be stronger than the present evidence against 
SMA. In this situation one could doubt that Borexino will 
be able to detect solar neutrinos at all. Therefore we also 
assume that 

2. SNO will see a distortion of the spectrum, and will 
contradict the SK spectral data.^ 

Assumptions 1 and 2 imply a very bad goodness-of-fit. A 
good SMA fit could be obtained if the SK collaboration 
would commit hara-kiri, admitting: 

3. serious faults in the SK solar neutrino results. 

In this case, we would be authorized to drop the the SK 
spectral data from the global fit. Under these three as- 
sumptions, we obtain the SMA future fit in fig. |l|. Borex- 
ino data select the SMA region over the other ones, but 
without favouring any particular corner there; this selec- 
tion is done by the other data, that prefer the region of 
large 9. We can consider a different possibility: if the CC 
rate measured at SNO is wrong (so that we drop it from 
the x^), but the SK spectral data are correct, the SMA 
range with small 9 would be selected. Incidentally, this 

t The first SNO spectral data indicate that this is not the case. 
The final SNO energy spectra are expected to be as significant as 
those of SK (SK will have more statistics, but the measurable recoil 
electron energy in SNO is more strongly correlated with the neutrino 
energy than in SK). In general, we have no reason to suspect that any 
solar neutrino data be wrong. However, if SMA were the solution of 
the solar neutrino anomaly, some of the present data that strongly 
disfavour SMA should be wrong. It is not true that Ax^ ~ 14 is 
less statistically significant than a direct \/T4 ~ 4 standard devia- 
tions evidence, because such a Ax'^ is obtained by merging several 
data or because oscillations have more than one free parameter. In 
fact, Borexino should measure a Beryllium rate 4 standard deviations 
lower than the one predicted by the current best-fit LMA solution, 
in order to make again SMA the 'best' fit solution. 
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Figure 4: Values of the Boron and Beryllium fluxes re- 
quired at la by the four measured solar neutrino rates in 
absence of oscillations, compared to the 68% C.L. ellipse 
of solar model predictions. 



tension between the preferred 9 range reflects once again 
the tension between the data, once we assume SMA. 

In any case, if something like this were to happen, we 
would certainly need new experiments to confirm it. The 
pp rate should be consistent with no oscillation, because 
SMA predicts Pee ~ 1 at very low energy. In order to make 
the solar neutrino anomaly credible, we would need to have 
a spectral measurement at low energy, aimed at revealing 
the sharp SMA transition from Pee « 1 to Pee ~ 0, for 
energies around the value in cq. (|^). 



and useful solar-model-indcpcndent information can be de- 
duced already from existing data |^ . 

The zero Beryllium flux required by no oscillations is 
disfavoured by helioscismology [|40[ , by simple physics con- 
siderations, and by recent results. Indeed, solar models |Q 
have been recently validated by two facts. The determi- 
nation of Sn (strength of the p^Be ^67 reaction) has 
been improved |4^ , and as a consequence the correla- 
tion between the flux of Boron and Beryllium neutrinos, 
visible in fig. ^ will tighten (both neutrino fluxes are con- 
trolled by the a^He — *■ '^Be7 reactions; thence the rela- 
tive S'34-factor remains the dominant, common source of 
uncertainty). Furthermore, the SNO/SK measurement of 
the Boron flux Q], is in agreement with the solar-model 
prediction, within the ^ 20 % error. These two facts sug- 
gest that the calculated Beryllium flux should not be too 
wrong^ Certainly, the solar model independent considera- 
tions will become cogent after the Borexino measurement 
of the Beryllium flux. 

Since the prediction of the pp flux is very accurate |Q 
^pp = (5.96 ± l%)10^°/cm^s, its determinations will also 
have an important impact on these analyses. Indeed, a 
well known simple (though, less accurate) argument leads 
to the same value of the flux. We assume that the so- 
lar energy comes from nuclear reactions that reach com- 
pletion, and that the sun is essentially static over the 
time employed by photons to reach the surface. The to- 
tal luminosity of the sun, Kq = 8.53 • 10" MeV cm'^ 
s^^ at the earth, determines its total neutrino luminos- 
ity as Kq = {Q/2- {E^J)^o.; Q = 26.73 MeV is 
the energy released in the reaction Ap + 2e ^ ^e -|- 2i/e, 
and the sum extends over the various components (a = 
pp,pepjBe,^^'N,^^0,^'^,^,hep). Neglecting {E^J and 
considering only the dominant pp flux, one obtains <E>pp w 
2Kq/Q = 6.4 • IQio/cm^s, that is only 7% off. Therefore, 
a measurement of the pp rate, even with modest accu- 
racy, will provide strong solar-modcl-indepcndent evidence 
against (or for) the extreme possibility of no oscillations. 



No oscillations 

The hypothesis of no oscillations can be reconciled with the 
data if the Boron and Beryllium fluxes are very different 
(smaller) from solar model predictions, as shown in fig. |^; 
see and ref.s therein. Even in this case, the rates mea- 
sured by the various experiments are not well compatible 
between them. But, as shown in fig. ^, there is a sig- 
nificant (accidental?), partial overlap of the experimental 
values.^ This is why solar- model independent considera- 
tions cannot strongly disfavour the no-oscillations hypoth- 
esis — whereas generically such accident does not happen 

tif GNO will reduce the error (or the central value) of the Gallium 
rate, the crossing in fig. U will happen at negative unphysical values of 
the Beryllium flux. The discrepancy between the SK and SNO bands 
is the well known solar-model-independent evidence for appearance 
of u^^T neutrinos. 



Sterile neutrinos 

The LSND anomaly can be used to argue for additional 
sterile neutrino(s); however, the evidence against such neu- 
trino(s) is now stronger than the evidence for them. For 
example, the '2-|-2' best solution found in where LSND 
is explained has a worse global fit than the 3-neutrino so- 
lution where it is not (and consequently than the best 
'3-1-1' solutions, where the sterile neutrino is of little use 
for LSND Q). Indeed, the relevant global is: 

2 2 2 2 

Xglobal = Xsun + Xatm + XlSND" 

§ Conversely, nuclear data and existing measurements |^ tend 
to suggest that 534 (and thence the Beryllium flux) is on the large 
side of the predicted value. 
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Using one sterile neutrino, the best fit gives pq | 

Xgiobai = 40 + 36 + = 76 with 9 parameters. 

Without using the sterile neutrino, the is 

Xgiobai = 37 + 29 + 10 = 76 with 5 parameters. 

CP-violating phases are not counted as relevant parame- 
ters. The Xgiobai of does not take into account data 
from SNO; from Karmen; from nucleosynthesis; from un- 
published atmospheric SK data about the total 7r° rate 
(and recent K2K cross section determinations) and about 
their zenith angle dependence . Each one of these data 
further disfavours the best-fit '2-1-2' solution of |4^, which 
has a large sterile component in atmospheric oscillations. 
From a theoretical point of view, '2-1-2' schemes need a very 
special arrangement of mixing angles. This pattern can be 
obtained from a 'pseudo-Dirac' mass matrix, which can be 
justified by a broken U(l) — Lr — Lg symmetry, and 
implies quasi-maximal mixing in atmospheric oscillations. 

If MiniBoone will confirm LSND, it will be interesting 
to consider solar oscillations into a mixed sterile/active 
neutrino. What would be the impact on the present anal- 
ysis? The near-future prospects concerning the determi- 
nation of the solar parameters 9 and Am^ are almost un- 
changed. This is an exact statement in the LMA region, 
since KamLAND cannot distinguish if reactor D,, disappear 
into active or sterile antineutrinos. The most relevant new 
issue is that solar oscillations would depend on a third 
parameter (other than Am^ and 0), that quantifies the 
sterile component in solar oscillations (for a precise defi- 
nition, see e.g. @). In this case, it would be interesting 
to supplement the NC/CC measurement of the Boron flux 
performed at SNO and SK with a NC/CC measurement 
of the pp flux, because it is accurately predicted by solar 
models. Furthermore, it could be convenient to obtain a 
CC measurement of the Beryllium flux, which could be 
combined with the result of Borexino (KamLAND) . 

3 Summary 

In the near future KamLAND or Borexino should identify 
the true solution of the solar neutrino problem, if it is 
due to oscillations. Depending on the actual value of the 
oscillation parameters Am^ and 9, the future situation will 
be very different, and will correspondingly require different 
new experiments. We summarize the various possibilities 
below (see the text for a more detailed discussion): 

• LMA. KamLAND or sub-KamLAND wiU measure 
Am^ with few per- cent accuracy. Even the mixing 
angle 9 can be determined reasonably well by reactor 
experiments; it will be a real challenge for sub-MeV 
experiment to improve on these measurements. 



• LOW. Borexino will see day/night effects and mea- 
sure Am^ with 10% accuracy. A measurement of the 
pp rate would be useful for determining 9 and a mea- 
surement of day/night effects in pp neutrinos could 
help in determining Am^. 

• LOW/QVO boundary. No unmistakable oscillation 
effect will be found, but all other solutions will be 
excluded. A measurement of day /night effects (that 
are largest for the lowest-energy pp neutrinos) would 
be crucial. 

• QVO and VO. Borexino will see seasonal effects and 
measure Am^ with few per-mille accuracy, that can 
be improved with more statistics. A measurement of 
the pp rate would be useful. 

We comment also on strongly disfavoured possibilities, that 
could strike back again, if near- future experiments will con- 
tradict some combination of established data: 

• SMA. Will become again the best solution if Borex- 
ino finds almost no solar neutrinos (because all Beryl- 
lium Ve get converted); one has to assume e.g. that 
the SK spectral data are not correct. A spectral mea- 
surement around E^, ~ MeV would provide a crucial 
signal. 

• No oscillations. Could become the 'best' solution if 
Borexino finds no solar neutrinos (because the Beryl- 
lium flux is much smaller than what solar models 
and helioseismology tell us — this possibility looks 
very remote) and if the CC/NC rate measured at 
SNO/SK is incorrect. The measurement of the pp 
rate would be of essential importance. 

• JustSo^. Will become the best solution if the CC/NC 
rate measured at SNO/SK will change (contradicting 
present data) indicating no oscillation, and if Borex- 
ino will flnd a somewhat suppressed flux of Beryllium 
neutrinos. A measurement of the pp spectrum would 
be crucial. 

Near-future experiments will allow us to deduce the pp flux 
from GNO data, without assuming oscillations or using 
solar model predictions, with < 15% uncertainty. 

Acknov^fledgments We thank E. Bellotti, C. Cattadori, 
N. Ferrari A. de Gouvca, A. lanni, M. Junker and S. 
Schonert for useful discussions. 

A A seasonal foi" Borexino 

The standard procedure employs a certain number A^bins 
of seasonal bins and a analogous to the one defined in 
eq. (0) for the day/night analysis. As is apparent from 
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fig. ^ the appropriate number of seasonal bins is propor- 
tional to the unknown value of Ato^. and many seasonal 
bins are necessary when Ato^ < 10^^ eV^ because A'osc is 
large. However, eq. (|^) cannot be applied if the number 
of events in each bin is not much larger than 1. In this 
situation, one should not employ Gaussian statistics; in- 
stead, the likelihood of the given measurement, given the 
oscillation parameters reads: 



n 

i=l 



:2_g-£i 



where the individual factors are the Poisson probabilities 
of having rii observed events in the i"* bin, in which 
events are expected. This function can be written as: 



TV„ 

= 2iVth - 2 ^ m ^ log + const, 

m ijn — l 

where iVth denotes the total number of expected events and 
the index im runs over the N^n bins with m = {1, 2, . . .} 
observed events, so that Nq + Ni + N2 + ■ ■ ■ ^ 

Poisson statistics simplifies in the limit of a very large 
number of seasonal bins (e.g. 1 bin per millisecond), so 
that ^ 1. In this situation there are only bins with or 
1 measured events, so that A'o + A^i = A'bins and 



2iV, 



th 



25:10: 



(6) 



If A'osc were large, this x^ may be useful to perform a sea- 
sonal analysis at Borexino, already after the first 6 months 
of data. 

We compare our x^ with other approaches. A stan- 
dard Fourier-transform analysis of Borexino seasonal data 
has been suggested in Equation (^) is, essentially, a 
non-standard type of transform, performed with respect to 
the event rate predicted by oscillations, rather than to the 
standard Fourier basis of 'sin' and 'cos' functions. This 
non standard choice minimizes the uncertainty on 9 and 
Am^. A similar idea has been proposed in (where a 
rough approximation to the spectrum of Beryllium neutri- 
nos has been employed) but its statistical meaning is not 
clear. Systematic and theorethical uncertainties can be 
easily taken into account, by writing the expected number 
of events as a function of 'nuisance' unknown parameters, 
in analogy with eq. (H). Our definition of the x^ is es- 
sentially the same as those employed in analyses of the 
SN 1987A data Wk- 



^The constant is irrelevant for parameter estimate. This does 
not provide a useful goodness-of-fit test, as any with too many 
bins. 



B Details of the computation 

Unfortunately, fitting solar neutrino data is a subtle issue: 
In order to perform inferences on the parameters of oscil- 
lation, one has to merge together many pieces of data and 
information (nuclear physics, solar models, matter effects 
in the sun and in the earth, various experiments, . . A We 
used information from many papers ||, ^, ^ ^, 31 , 1^, 



H, m, H, H, H, |0[ |l|, g ^. It is briefer to list here 
what is not included in our fit. Most of the other global 
fits have similar shortcomings. 

The spectrum of recoil electrons in SK and SNO is 
computed in the simplest approximation (e.g. neither one- 
loop effects nor photon emission (54j are included). One- 
loop corrections to the MSW effect are neglected. Seasonal 
Gallex/GNO and SAGE data are not included. The treat- 
ment of solar matter effects in the QVO region is not as 
precise as in |]3^ . When computing confidence levels, we 

approximate 1'-^ with a Gaussian function of Arn^ and 
B: correct frequentistic and Bayesian analyses po|| do not 
give a significantly different result. All above issues do not 
have significant effects. Alternative possible definitions of 
the x^ give A^^ values similar to those quoted in table 1. 
Note that the first digit of the Ax^ is significant. 

Our fit of SK data is based on table III of II . It only al- 
lows us to reproduce the total rate and the total day/night 
asymmetry (quoted in |^) with ~ \a accuracy. This small 
discrepancy is presumably due to the fact that the used 
by the SK collaboration includes a proper treatment of 
data about the background, so that the total rate is not 
the sum of the rates in each energy bin: the bins with 
higher energy are relatively more important, since they 
have less background. 
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